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Abstract: We stress the potential usefulness of renormalization group invariants. Espe¬ 
cially particular combinations thereof could for instance be used as probes into patterns of 
supersymmetry breaking in the MSSM at inaccessibly high energies. We search for these 
renormalization group invariants in two systematic ways: on the one hand by making use 
of symmetry arguments and on the other by means of a completely automated exhaustive 
search through a large class of candidate invariants. At the one-loop level, we find all 
known invariants for the MSSM and in fact several more, and extend our results to the 
more constrained pMSSM and dMSSM, leading to even more invariants. Extending our 
search to the two-loop level we find that the number of invariants is considerably reduced. 
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1 Introduction 

It has long been known that the Standard Model cannot be the final theory of particle 
physics. Issues such as the hierarchy problem or the absence of a satisfactory description 
of dark matter or gravity lead theoretical physicists to develop more fundamental theories. 
Supersymmetry is one of those theories, and the Minimal Supersymmetric Standard Model 
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(MSSM) will be one of the main focuses for the second run of the Large Hadron Collider. 
If signs of supersymmetry are indeed found, theorists will face the issue of figuring out 
the exact underlying theoretical description. In case of the MSSM, this problem includes 
finding the exact mechanism of supersymmetry breaking. If supersymmetry is a symme¬ 
try of nature, every particle must have the same mass as its superpartner. Since this is 
clearly not the case, supersymmetry must be broken. The MSSM accounts for several 
breaking mechanisms by incorporating all possible soft supersymmetry breaking terms in 
its Lagrangian. If any of the several possible breaking mechanisms is realized in nature, 
this will be signified by a characteristic unihcation of some soft supersymmetry breaking 
parameters at a high energy scale. 

These unihcations can be studied through renormalization group (RG) techniques. 
Typical strategies include evolving the values of measured parameters upward to the uni¬ 
fication scale (bottom-up), or choosing values of the parameters at the unihcation scale 
that are evolved downward (top-down). We will discuss a third method that makes use of 
RG invariant combinations of RG equations, called RG invariants. Using RG invariants to 
probe high-scale physics, has several distinct advantages over the bottom-up and top-down 
methods [1-3]. Finding these invariants can be very difficult though. In this paper we 
will discuss two different methods to hnd RG invariants. One method relates invariants to 
symmetries of the underlying theory and the other is based on computer algebraic tech¬ 
niques. Both methods will be applied to the MSSM, the diagonal MSSM (dMSSM), and 
the phenomenological MSSM (pMSSM) to hnd invariants up to two-loop order. 

The paper is organized as follows. In section 2 we describe the method of RG invariants 
and its merits. In section 3 we investigate whether some of these RG invariants can be 
related to symmetries of the underlying theory. For this underlying theory we will consider 
the MSSM, dMSSM, and pMSSM. In section 4 a generic computer algebraic method for 
hnding RG invariants is introduced. This method is used to hnd previously unknown 
invariants at one- and two-loop level within the aforementioned supersymmetric theories. 

2 Probing high-scale physics with RG invariants 

Let us consider a renormalized theory with a running parameter p{^). We dehne the 
corresponding /3-function as follows: 

/3(p) = (2.1) 

where t = logiQ{fj,/po). The energy scale /i is normalized by an arbitrary reference scale /tq 
to make the logarithm dimensionless. For quantum held theories, the one-loop /3-functions 
are polynomials of the parameters of the theory with rational coefficients. Higher order 
contributions contain additional factors of l/(167r^). 

How could RG invariants be used to probe high-scale physics? We will explain this 
through a toy system of one-loop /3-functions for the parameters v,w,x,y,z that closely 
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resembles some of the one-loop MSSM /3-functions; we define 


II 

CO 

(2.2a) 

f3{w) = {5w + 6x — 4y), 

(2.2b) 

I3{x) = {—w — 2x + Az ), 

(2.2c) 

/3(y) = v‘^{x + y-hz), 

(2.2d) 

(3{z) = {w — 2y + Qz). 

(2.2e) 


This system of /3-functions was built such that the parameter v resembles a gauge coupling 
parameter, while w, x, y, z resemble MSSM scalar masses. An RG invariant is an algebraic 
combination of parameters I, such that 

(2.3) 

As it turns out, two independent RG invariants can be constructed from the /3-functions 
given by eqs. (2.2a) - (2.2e); we define 

Ii = w + 3x — 2z, l 2 ^ X + 2y + z. (2.4) 

Supersymmetry breaking mechanisms typically predict the unification of scalar masses at 
some large, experimentally inaccessible energy scale. We can use the invariants in (2.4) to 
check whether the unification of w, x, y, z is realized in nature. Suppose that these scalar 
masses unify to the value s at some scale, then we would have 

h = 2s, h = 4s, (2.5) 

from which it follows that 

2h - /2 = 0 . ( 2 . 6 ) 

Such a relation between RG invariants is called a sum rule. Since Ii and I 2 are invariant 
under RG flow, sum rules such as eq. (2.6) remain true at all energy scales, if unification 
occurs. In particular, sum rules can be checked at the collider scale, potentially falsifying 
the assumed unification. 

The approach using RG invariants avoids a number of issues that the top-down and 
bottom-up methods suffer from. The top-down method requires knowledge of both the 
unification scale and value in order to evolve the parameters down to experimentally ac¬ 
cessible scales. Since theory does not predict these values with much accuracy, scans are 
typically performed over a range of scales and values. These scans can be computationally 
very time-consuming (depending on the number of unifying parameters), and are often not 
feasible. The use of RG invariants requires no knowledge of the unification scale or value 
whatsoever, so that this problem is avoided. 

The bottom-up method does not require knowledge of unifying scales or values either, 
but suffers from a different problem resulting from the numerical evolution of experimental 
input values. When parameters are evolved up to higher scales, any experimental errors 
are typically greatly enhanced. Figure 1 shows the running effects of a slight change of 
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one of the low-scale values of the parameters v,'w,x,y,z when w,x,y,z unify at some high 
scale. After the slight change, the evolution of the parameters no longer shows any sign of 
unification whatsoever. The use of RG invariants circumvents such numerical problems. 



Figure 1. Left: The evolution of the parameters v, w, cc, y, z according to their one-loop ,d-functions 
(2.2a) - (2.2e). The parameters w,x,y,z unify at t = 16 to the value 10. Right: The initial value 
of y at t = 2 is raised by 1%. 

Finally, it is usually not necessary to know all parameters of the theory when using RG- 
invariants as probes of high-scale physics. The /3-functions for most theories are heavily 
coupled differential equations, and hence evolving the relevant parameters of a theory 
usually requires the evolution of all other parameters of the theory as well. In our example 
the sum rule of eq. (2.6) does not involve the parameter v, which does not participate in 
the unification anyway. Both the top-down and bottom-up methods also require a value 
for V to numerically evolve the other parameters. 

Constructing RG invariants is generally a very non-trivial matter. To find invariants, 
we will not use the “standard” algebraical techniques that are used in for example [3]. 
Instead, we will consider two different, more efficient methods: one relies on symmetries 
of the underlying theory, and the other is based on computer algebraic techniques. Both 
methods will be applied to find invariants for the MSSM, the dMSSM, and the pMSSM. 

3 RG invariants from symmetries 

First, we will look at the construction of RG invariants from a perspective that involves 
symmetries of the underlying theory. The existence of such a relation between symmetries 
and invariants was suggested by [2] in the context of the pMSSM. As yet, no compelling 
arguments or proof of this has been presented though. First we will consider the MSSM, 
then we will study a “flavor-diagonal” version of the MSSM (the dMSSM), and, finally, 
the pMSSM will be discussed. All nomenclature regarding MSSM, dMSSM, and pMSSM 
fields and parameters is defined in appendices A, B, and C respectively. 
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3.1 The MSSM 


In this subsection we attempt to construct RG invariants in the MSSM. So far, RG in¬ 
variants have only been constructed for heavily constrained supersymmetric models, such 
as the pMSSM [2, 3]. We will attempt to construct invariants for the full 168-parameter 
MSSM. 

3.1.1 Known invariants 


We expect that some invariants that have been found in the pMSSM also exist in the 
MSSM, since the simplifications of the pMSSM with respect to the MSSM mainly apply to 
the family sector, and, for example, not at all to the gauge sector. The following invariants 
that have been found in the pMSSM are also invariants in the MSSM, which can be checked 
easily with the /3-functions in appendix A.4. In the combined gauge and gaugino sectors, 
we can construct three invariants:^ 


_ Ml _ M2 _ M3 

9 ^ 9^ 9i 


and in the pure gauge sector we have 

1 


h = 


/2 

9 9 


11 
w ■ 


7-3 11 

-'s = ^ + W 
9 9i 


(3.1) 


(3.2) 


Finally, the invariant that involves the quantity S turns out to be an invariant of the MSSM 
as well: 

(3.3) 

9 

Thus, six invariants that we already know from the pMSSM carry over trivially to the 
MSSM. 

3.1.2 New invariants 

As we know from [2, 3], all remaining invariants in the pMSSM (there exist eight more be¬ 
sides Ii,..., Ig) involve scalar masses only, or combinations of scalar and gaugino masses.^ 
For this reason, we will now focus on these sectors in the MSSM to see if more invariant 
quantities can be constructed. 

Let us now try to construct invariants from the one-loop /3-functions of the soft scalar 
masses and gaugino masses that are listed in appendix A. 4. The first thing we note when 
we look at the /3-functions for the soft scalar masses, is that no invariants can possibly be 
constructed for the off-diagonal terms of the sfermion mass matrices. The reason for this is 
that the order in which the family space matrices (i.e. the Yukawa and trilinear coupling 
matrices) appear is different for all five /3-functions of the sfermion mass matrices. In other 
words, there are simply too many different structures present. In order to avoid this, we 

^Strictly speaking there are six invariants: these three, pins their complex conjngates. 

^Actually, in the pMSSM two more invariants exist that involve the Higgs mixing parameters fi and b, 
but as argued in [3] these are useless to probe high-scale physics models. For completeness, though, they 
are listed in appendix D. 
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need to work with terms that are insensitive to this order of matrices. Hence, a logical step 
would be to consider the /3-functions for the traces of the sfermion mass matrices. These 
/3-functions are given by 




Tr(m|) 


/3 



/3 [TV(m|)] 





/3 [Tr(m|)] 






- - 18ff2|M2|2 - 32g^,\Ms\^ + 

+ 2Ti (mjjy^yu + m^^ylyu + ylmlyu + ata«) 
+ 2 TV (mjj^ylyd + m^^ylyd + y\m^^yd + a^aa^ , 

- 24y|^c/'2|Mi|2 - ISg^lMap + GYij'^S 

+ 2Ti {m\jjjlye + rn^ylye + ylm-ye -V a^ae) , 

- 2AYi.J^\Mi\^ - 3252|M3|2 + 6Y^*J^S 

+ 4TV (mjj^yly^ + m^^ylyu + ylmlyu + , 

- 24Tj, - 32<72|M3|2 + OTj. g'^S 

+ 4 TV (mjj^ylyd + m^-y^yd + y^m^yd + a^a^) , 

- 24y4ff'2|Mi|2 + 6Ye*J^S 

-h 4 TV (mjj^ylye + m^yly^ -V ylm^y^ -V a^ae) , 

- 3Y^J^\Mi\^ - 6ff2|M2|2 + 2YhJ^S 

+ 6TV {m]j^ylyu + m^~ylyu + ylmlyu + a^a^) , 

- + 2YhJ^S 

+ 6 TV {m]j^ylyd + m?~ylyd + ylm^yd + a^a^) 

-F 2 TV (mjj^ylye + m^yly^ + ylm^y^ -V , 


(3.4a) 


(3.4b) 


(3.4c) 


(3.4d) 


(3.4e) 


(3.4f) 


(3.4g) 


where, for convenience, we have also added the /3-functions for and (which 

already involved traces). 

In the MSSM, all interactions in family space are described by the Yukawa terms in 
the superpotential and the soft supersymmetry breaking trilinear terms (cf. appendix A. 2). 
In case that the order of the family space matrices does not matter by taking a trace, then 
these interactions give rise to three different trace structures in (3.4a) - (3.4g), labeled by 
u, d, and e. These trace structures each belong to one of the three Yukawa and trilinear 
interaction terms in the superpotential and the soft breaking Lagrangian respectively, and 
each of those terms involves a unique set of three scalar fields. For example, the terms in 
the MSSM Lagrangian that give rise to the trace structure that features the label u are 


ATFmssm = u\^yu{QLT{Hu)a, A.iftrii. = -u^RCLu{QLT{Hu)a + h.c. (3.5) 


Note that the coefficients of the three trace structures are different for the various 
scalar masses. This is because in the /3-functions no traces have been carried out over the 
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gauge degrees of freedom of the corresponding scalar fields. How is this to be understood? 
Let us consider the trilinear interaction between the fields ur, Ql, and as well as the 
one-loop corrections to the scalar propagators that this interaction gives rise to. These 
three fields can all occur as external fields or inside loops. The gauge degrees of freedom 
that are “closed” inside the loops (i.e. those gauge degrees of freedom that the external 
fields do not possess) are summed over. Say we take ur to be the external held, then 
there is an SU(2) doublet degree of freedom inside the loop that has to be “traced over”, 
giving a factor of 2. If we had taken Hu to be the external held, then a trace over SU(3) 
degrees of freedom inside the loop would have resulted, giving a factor of 3, etc. Hence, if 
we multiply the /3-functions (3.4a) - (3.4g) by the factors that result from summing over 
the gauge degrees of freedom of the external scalar helds, then all three trace structures 
each get exactly the same coefficients. 

How many RG invariants do we expect to hnd in the soft scalar and gaugino sectors? 
We have ten equations (seven /3-functions for the scalar masses and three for the gaugino 
masses) that contain seven different structures (three trace structures, three gaugino masses 
and S). Having ten equations to eliminate seven different terms should give 10 — 7 = 3 
independent RG invariants. 

To construct invariant quantities from these /3-functions, we can hrst try to get rid 
of the three different trace structures. To cancel these structures, we could consider a 
linear sum of the /3-functions (3.4a) - (3.4g), appropriately multiplied by numbers of gauge 
degrees of freedom, such that the coefficients of this sum add up to zero for each trace 
structure. Now we can use the fact that each trace structure corresponds to a unique 
combination of three scalar fields in the Lagrangian. For the cancellation to take place, we 
should assign quantum numbers to these groups of three fields that each add up to zero, 
which is equivalent to saying that the quantum number should be conserved by all family 
space interactions. This means that to cancel the Yukawa and trilinear contributions to 
the /3-functions, we should consider U(l) symmetries of the MSSM family sector. More 
specifically, for any quantum number Q that pertains to a symmetry U(1 )q of the MSSM 
family sector, the /3-function of the quantity 



with the sum running over all scalar fields 0, does no longer contain the three trace struc¬ 
tures. Note that sums like the one in (3.6) over fields that occur inside family space traces 
are implicitly understood to be over all gauge degrees of freedom of the gauge multiplets.^ 
What could the quantum number Q be? As is discussed in appendix A. 3, the relevant 
U(l) symmetries correspond to the quantum numbers weak hypercharge (Y), baryon num¬ 
ber (B), lepton number (L), and X. Table 1 provides the quantum numbers Y,B,L,X for 
all scalar helds in the MSSM. The presence of exactly/owr independent U(l) symmetries in 
the family sector of the MSSM can be explained as follows: from table 1 we infer that this 

®For the field Ql, for example, such a sum would give a factor of 2 for SU(2) and a factor of 3 for SU(3) 
degrees of freedom, yielding a total multiplication by 6. 
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sector of the theory is constituted by seven scalar multiplets. We can regard each set of 
quantum numbers pertaining to a given symmetry as being a vector in a seven-dimensional 
vector space. This vector is subject to three independent symmetry constraints coming 
from the interaction terms in the superpotential and the soft trilinear terms (the terms in 
eq. (3.5) give one of these constraints). This means that we can construct 7 — 3 = 4 linearly 
independent vectors in this space, i.e. four sets of quantum numbers each corresponding 
to a different U(l) symmetry. 


Spin 0 

Y 

B 

L 

X 

Ql 

1 

6 

1 

3 

0 

1 

Ll 

1 

2 

0 

1 

1 

y*R 

2 

3 

1 

3 

0 

1 

^*R 

1 

3 

1 

3 

0 

1 


1 

0 

-1 

1 

Hu 

1 

2 

0 

0 

-2 

Hd 

1 

2 

0 

0 

-2 


Table 1. The quantum numbers Y,B,L,X for all MSSM scalar fields. 


We now know how to get rid of the three family space trace structures in eqs. (3.4a) 
- (3.4g), but what about the other terms present? Can quantum numbers also be used to 
eliminate the remaining structures (i.e. S and the three absolute squared gaugino masses 
|Mip, IM 2 P, IM 3 P)? For any Q, we have 


/3 


Tr 









(3.7) 


where d denotes the scalar weak isospin doublets and t the scalar color triplets. Note that 
sums without any family space traces involved are implicitly understood to be over all 
families as well (besides the gauge degrees of freedom).'^ From eq. (3.7) it directly follows 
that to eliminate S, |Mip, IM 2 P, IM 3 P respectively, we must have: 


= 0, = J]Qd = 0, ^Qt = 0. (3.8) 

(f) <p d t 


The latter three sums over charges are reminiscent of mixed anomaly cancellations of the 
charge Q with the separate gauge groups that we know from for example the Standard 
Model. 

^This simply amounts to an additional factor of 3 for all sfermions. 
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Which quantum numbers actually satisfy the requirements in (3.8)? This is summa¬ 
rized in table 2. From this table it follows that to eliminate S, suitable quantum numbers 
would be 3B + L, 11{B — L) — 8Y, and X. For the cancellation of |Mip, we could use Y, 
B — L, ov 16B + 3X, while for IM 2 P the quantum numbers Y, B — L, and 8 B — 3X are 
suitable. For [Map to cancel, we could pick Y, B, or L. Naturally, linear combinations of 
these quantum numbers also work. 


Quantity 

Sum 

Y 

B 

L 

X 

S 


11 

2 

-6 

0 

|Mi|2 


0 

3 

2 

3 

2 

8 

|M2|2 

Ed Qd 

0 

6 

6 

16 

CO 

to 

Et Qt 

0 

0 

0 

36 


Table 2. The evaluation of the sums that are related to the elimination requirements (3.8) for the 
quantities S, |Mi^ 2 , 3 p, given for the quantum numbers Y,B,L,X. The outcome 0 for a certain 
quantum number indicates that this quantum number is suitable for eliminating the corresponding 
quantity from the /3-frmctions (3.4a) - (3.4g). 


From eq. (3.7), as well as from the /3-functions for the gaugino masses (A.20) and 
S (A.23), it follows that any one-loop RG invariant I in the MSSM that involves scalar 
masses and gaugino masses, is of the following form: 

iq = tv(y1 - n En E^ Qdmf 

^ <p (p (p d 

(3.9) 

t 

where Q is any quantum number that is preserved by all MSSM family space interactions 
(i.e. Q must be a linear combination of Y, B, L, X). 

Now we are ready to construct RG invariants using table 2 and eq. (3.9). From table 
2 we infer that the quantum number ll(i3 — L) — 8 Y cancels S and all gaugino masses, 
as this linear combination of Y, B, L vanishes for each row in the table. This leads us to 
define the following RG invariant: 

77 = Tr^^ {IIB^ - IIL^ - 8 Y^) 

^ 0 

= Tr ^14m~ — 14m~ -|- 5m~ — 19m3m|^ — 8 na%^ + (3.10) 

If we wish to construct an invariant where both S and | M 3 1 ^ are eliminated by a symmetry. 
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then we could use the quantum number 3B + L. We define 

Is = Tr ( ^ 

^ <p ^4) d 

= Tr (^6m~ + 2m~ — 3m~ — 3m?~— m|^ — ^|Mip + 36|M2p. (3-11) 

For the third and last independent invariant in this sector, let us pick the quantum number 
X (which only cancels S) and define^ 

/, s Tt( X^ml) + A ^ v?^ E ^ E 

^ 0 0 d t 

= Tr (^6m~ + 2m| + 3m~ + 3m~+ m|^ - - 4m|^^ + ^|Mip 

+ 24|M2p-32|M3|2. (3.12) 

3.2 The dMSSM 

Let us consider a constrained version of the MSSM, the so-called dMSSM, where all matrices 
in family space (i.e. the sfermion mass matrices and the Yukawa and trilinear coupling ma¬ 
trices) are taken diagonal (cf. appendix B for a more extensive discussion of the dMSSM), 
and let us again focus on the soft scalar and gaugino sectors to find RG invariants. In this 
particular model, we have twenty equations (fifteen /3-functions for the sfermion masses, 
two for the Higgs masses, and three for the gaugino masses) containing thirteen different 
structures (nine structures coming from the diagonal components of the Yukawa and tri¬ 
linear coupling matrices, three gaugino masses, and S). Eliminating only thirteen different 
terms using twenty equations would result in 20 — 13 = 7 RG invariant quantities for this 
simplified model, on top of R,..., Jg. 

The “traced” /3-functions for this simplified model are again given by eqs. (3.4a) - 
(3.4g) and, as we discussed in the previous subsection, give rise to three independent RG 
invariants of the form (3.9). Thus, the invariants Ij, /s, Lg are also invariants in this model. 
How can we construct the remaining four invariants? Can we again benefit from symmetry 
arguments? 

As we have taken all family space matrices diagonal, there is no longer any flavor mixing 
present. In other words, the three (s)fermionic generations have completely decoupled. This 
means that the baryon and lepton numbers are separately conserved for each generation, i.e. 
B and L can now be split up into Bi,B 2 , B^ and Li, L 2 , Ls respectively.® Table 3 provides 
all quantum numbers that pertain to (independent) U(l) symmetries of the dMSSM. The 
completeness of this list can be shown in the same way as we did for the MSSM: from table 
3 and appendix A.2 we infer that we have seventeen gauge multiplets that are subject to 

®Even though we have four symmetries at hand, only three independent RG invariants can be con¬ 
structed. This is because taking Q oc Y gives Iq = 0. In fact, the quantum number Y has already been 
used for the invariant that involves the quantity S (cf. eq. (3.3)). 

®Note that two linear combinations of these six quantum numbers are equivalent to B and L, namely 
the combinations Bi + B 2 + Bs and Li + L 2 + L 3 respectively. 
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three constraints per generation. Hence there must be 17 — 3 x 3 = 8 independent U(l) 
symmetries in this model. 


Spin 0 

Y 

Si 

S 2 

S 3 

Si 

S 2 

S 3 

X 

{Ql)i 

1 

6 

1 

3 

0 

0 

0 

0 

0 

1 

{Ql)2 

1 

6 

0 

1 

3 

0 

0 

0 

0 

1 

{Ql)3 

1 

6 

0 

0 

1 

3 

0 

0 

0 

1 

(^l)i 

1 

2 

0 

0 

0 

1 

0 

0 

1 

(Ll)2 

1 

2 

0 

0 

0 

0 

1 

0 

1 

(Ll)3 

1 

2 

0 

0 

0 

0 

0 

1 

1 


2 

3 

1 

3 

0 

0 

0 

0 

0 

1 


2 

3 

0 

1 

3 

0 

0 

0 

0 

1 


2 

3 

0 

0 

1 

3 

0 

0 

0 

1 


1 

3 

1 

3 

0 

0 

0 

0 

0 

1 

(4)2 

1 

3 

0 

1 

3 

0 

0 

0 

0 

1 

{d*Rh 

1 

3 

0 

0 

1 

3 

0 

0 

0 

1 

(^k)i 

1 

0 

0 

0 

-1 

0 

0 

1 

(6^)2 

1 

0 

0 

0 

0 

-1 

0 

1 

{^Rh 

1 

0 

0 

0 

0 

0 

-1 

1 

Hu 

1 

2 

0 

0 

0 

0 

0 

0 

-2 

Hd 

1 

2 

0 

0 

0 

0 

0 

0 

-2 


Table 3. All quantum numbers that pertain to U(l) symmetries of the dMSSM family sector. The 
dashed lines separate the different family multiplets. 


Since there are four new independent symmetries in the dMSSM with respect to the 
MSSM, four additional invariants of the form 

iq = Y1 - n E ^ E ^ E ^ E 

(p (f) (f) d t 

(3.13) 

can be constructed. Clearly, the quantum numbers Bi — B 2 , Bi — B^, Li— L 2 , and Li — L 3 
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automatically eliminate both S and all gaugino masses. We define 


^10 = + m|, (3.14a) 

<t> 

Ill = = 2m|^ - - m| - 2m|^ + mig + m|, (3.14b) 

<i> 

h 2 = = 2m|^ - m|^ - 2m|^ + (3.14c) 

4 > 

/i3 = “ ^30)m^ = 2m|^ - m2 - 2m|^ + m^. (3.14d) 

<i> 


3.3 The pMSSM 

For the pMSSM, with respect to the dMSSM, there are a couple of additional constraints: 
the first two generations of sfermions are mass degenerate, and the Yukawa and trilinear 
coupling matrices only have non-zero entries for the third generation sfermions (cf. ap¬ 
pendix C for a more extensive discussion of the pMSSM). This means that for the pMSSM 
we have fifteen equations (ten /3-functions for the sfermion masses, two for the Higgs masses, 
and three for the gaugino masses) to eliminate seven different structures (three structures 
coming from the (33)-components of the Yukawa and trilinear coupling matrices, three 
gaugino masses, and S), which should yield 15 — 7 = 8 invariant quantities on top of 
Ii,... ,Iq. Indeed, for the pMSSM eight RG invariant quantities have been constructed in 
[2, 3] that involve scalar and gaugino masses only. Can we construct these invariants also 
based on our approach involving symmetries? 

Again, three independent invariants can be constructed that are of the form (3.9), thus 
It, I%,Ig trivially carry over to the pMSSM. Of the invariants Iio,..., /is, only In and Rs 
also exist in the pMSSM (the invariants Iio and I 12 vanish due to the mass degeneracy of 
the first two sfermionic generations).' 

Are there additional symmetries in the pMSSM with respect to the flavor diagonal 
model of the previous subsection? The first thing to note is that the first and second 
generations of sfermions are completely identical in the pMSSM, which means that we 
need to consider twelve gauge multiplets only (ten sfermionic multiplets and two Higgses). 
As the Yukawa and trilinear coupling matrices only have non-zero entries for the third 
generation sfermions, no first and second generation sfermions feature in any of the family 
space interactions. This in turn means that any set of quantum numbers for the first two 
generations automatically corresponds to a symmetry of the pMSSM. The only constraints 
on allowed sets of quantum numbers (corresponding to symmetries) arise from the third 
generation sfermions. Similar to the MSSM, we have three constraints coming from the 
interaction terms in the superpotential and the soft trilinear couplings. This means that 
there exist 12 — 3 = 9 independent U(l) symmetries in the pMSSM family sector. Table 4 
provides the quantum numbers that correspond to nine independent U(l) symmetries of 

^Another way to look at this is that in the pMSSM the quantum numbers Bi , B 2 and Li, L 2 are each 
equivalent to each other, which means that there are only two new symmetries (and thus two new invariants) 
with respect to the MSSM. 
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the pMSSM. 


Spin 0 

Yi 

Y 

Bi 

Bs 

Li 

Ls 


Xlg 

X 

{Ql)i 

1 

6 

1 

6 

1 

3 

0 

0 

0 

0 

1 

1 

{Ql)3 

0 

1 

6 

0 

1 

3 

0 

0 

0 

0 

1 

(Xl)i 

1 

2 

1 

2 

0 

0 

1 

0 

1 

0 

1 

(Xl)3 

0 

1 

2 

0 

0 

0 

1 

0 

0 

1 


2 

3 

2 

3 

1 

3 

0 

0 

0 

0 

1 

1 


0 

2 

3 

0 

1 

3 

0 

0 

0 

0 

1 


1 

3 

1 

3 

1 

3 

0 

0 

0 

0 

1 

1 

('^r)3 

0 

1 

3 

0 

1 

3 

0 

0 

0 

0 

1 


1 

1 

0 

0 

-1 

0 

1 

0 

1 

(e«)3 

0 

1 

0 

0 

0 

-1 

0 

0 

1 

Hu 

0 

1 

2 

0 

0 

0 

0 

0 

0 

-2 

Hd 

0 

1 

2 

0 

0 

0 

0 

0 

0 

-2 


Table 4. All quantum numbers that pertain to U(l) symmetries of the pMSSM family sector. The 
dashed lines separate the different family multiplets. 


Let us now construct the three remaining invariants in the pMSSM, making use of the 
newly available symmetries and eq. (3.13). A suitable quantum number to cancel both S 
and all gaugino masses is 10(ili — Li) — 8Yi. For the other two, let us pick Xi^ and Xig. 
We define 

/i4 = — lOLi,^ — 8Yi0)m^ 

0 

= 12m~ — 12m| + 6m~ — I8m~ +2m~. (3.15a) 

/i5 = ^ ^ ^ XuM 

(p (j) d 

- 2m|^ + + d-Aff + 3 M 2 ) (3.15b) 

ho ^ E + n ^ 5 E xmM - ^ E 

(j) (p d t 

= 6m~ + 3m| + 3m~ + + 9M| - ^M|. (3.15c) 

^ “1 3 3 

The invariants found for the pMSSM are consistent with the ones derived in [2, 3]. 
®The quantum numbers in this table, of course, are not unique. 
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4 RG invariants from computer algebraic techniques 

Next, we approach the problem of finding RG invariants from a computer algebraic angle. 
The goal of this approach is to find a method that can in principle be applied to any set of 
/1-functions, for any theory. By letting a computer do the heavy lifting, we will not have to 
rely on any properties of the underlying theory as in the case of the previous method, but 
we will instead be limited by the available computational power. To develop the method, 
we consider two specific forms of one-loop invariants before extending the method to higher 
loop orders. 


4.1 Monomial invariants 

Let us first consider the simple class of monomial invariants. A monomial invariant M has 
the following form: 




(4.1) 


2=1 


where Xi (with z = 1,..., n) are running parameters and a G Z”. The requirement for RG 
invariance simply reads 


^ = 0 , 


dt 


(4.2) 


2 = 1 


from which it follows that for all values of the parameters Xi, we must have 


A ai/3{xi) ^ ^ 


(4.3) 


2 = 1 


To see how this works in practice, let us consider a simple toy system for two parameters 
X and y with 

f5{x) = xy+ 3xy‘^, P{y) = 2y‘^ + 6y^. (4.4) 

From requirement (4.3), it follows that 

ax/3{x) ay/3{y) _ , o ^ i ^ 2 _ n 

- 1 - — [cix + 2aj^) y -|- (Sox -|- Qdy) y — 0, (4.5) 

X y 


or in matrix form: 



(4.6) 


Finding a RG invariant has now reduced to solving for the nullspace of a linear system of 
equations. We find that a = (2, —1) spans the nullspace, which leads to the invariant 


I = 


y ' 


(4.7) 


Note that the existence of only a single invariant already implies the existence of an 
infinite amount of invariants, since a product of invariants is also an invariant. However, 
each of these invariants must solve eq. (4.3), and must therefore be included in the nullspace 
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of the linear system of equations. By finding a basis vector for this nullspace, we are 
effectively including all of these solutions. A higher-dimensional nullspace would imply 
the existence of multiple independent invariants. Any products of these invariants are also 
invariant, but their existence is again implied by the linearity of the problem. 

4.2 Polynomial invariants 

Let us now consider polynomial invariants. Such an invariant P has the following general 
form: 

m n 

i=i *=i 

The powers Ojj of the parameters Xi now in fact form a matrix a G We have also 

introduced C G Z™' as a vector that contains the coefficients for the separate monomial 
terms. The invariance requirement amounts to 



Unlike the case for monomials, we cannot factorize the invariant itself and we are left 
with a highly nonlinear equation in both the unknowns a and C. In addition to the fact 
that considering products of invariants yields an infinite set of solutions to eq. (4.9), now 
also linear combinations contribute to this issue. Clearly, a method for finding polynomial 
invariants must be able to deal with both of these sources for ending up with an infinite 
number of solutions. 

To fix the issue for products of invariants, we introduce the concept of dimensionality. 
Let us consider the following toy system for two parameters x and y, with 

(3{x) = 2xy + Wy^, (3{y) =-x - y"^. (4.10) 

We now assign a dimension to the parameters of this system. If we set dim(x) = 2 and 
dim(y) = 1, it follows that dim(/3(x)) = 3 and dim(/3(y)) = 2. Thus, we find that for all 
parameters xt, dim(/3(a:i)) — dim(xi) = c, where c is a constant. As a consequence, for any 
monomial M, 

= dim(M)-|-c. (4.11) 

In particular, if we consider two monomials M and M' with dim(M) ^ dim(M'), then eq. 
(4.11) implies that d\m{dM/dt) ^ diva.{dM'/dt). As a consequence, the monomial terms 
in dM/dt must be different from those in dM'/dt. Hence, if both M and M' are to be 
included in an invariant, they must be part of two separate groups of monomials that are 
separately invariant. But then we are considering a linear combination of invariants, which 
we are trying to avoid. Therefore, whatever the rest of the method is, no results are lost 
by considering invariants consisting of monomials of the same dimensionality. In addition, 
this deals with the issue of having an infinite set of solutions to eq. (4.9) due to considering 
products of invariants, since the product of two dimensionful invariants has a dimensionality 
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that is different from the two original invariants. Since multiple dimensionalities can be 
assigned to the parameters of a theory, the dimensionality of a monomial is in general a 
vector d. For a system of /3-functions with r distinct dimensionalities (i.e. d G Z^), we 
define a set of monomials of the same dimensionality as follows: 


Mp{d) 



a G Z”, G {1,..., r} : dim^ 



= di, 0 < tti < p 


(4.12) 


The restriction 0 < a* < p is simply included to ensure that Mp{d) is a finite set.® 
For the toy system of x and y, for example, we have 




(4.13) 


A candidate invariant can now be built by using the monomials in this set: 

S 

PliC) = ^C,Mi, (4.14) 


where Mj G Aip{d), and s denotes the size of Aip{d). Requirement (4.9) now gives 


0 = Ci2x/3{x) + C2y‘^l3{x) + C'22xy/3(y) + CsAy^Piy) 


= (IOC 2 - 4 C 3 ) y^ + ( 2 OC 1 - 4 C 3 ) + (4Ci - 2 C 2 ) x^y. (4.15) 


As this equality must hold for all values of the parameters x and y, we again recognize a 
linear system of equations: 


/ 0 10 -4 \ 
20 0 -4 

^4-20) 


[cA 

C 2 = 0. 

\Cs) 


(4.16) 


We find that the vector C 
invariant 


(1,2,5) spans the nullspace of this matrix, leading to the 


I = x^ + 2xy‘^ + 5y^- 


(4.17) 


Note that the reduction to a linear system of equations in the coefficients Cj auto¬ 
matically takes care of the problem of having an infinite amount of solutions to eq. (4.9) 
due to linear combinations of invariants. Those linear combinations are in fact just linear 
combinations of the vectors C. By finding the nullspace for a system of equations like 
(4.16), we handily deal with all problems of eq. (4.9) while maintaining a simple procedure 
that can be performed by a computer with ease. This method has been implemented in a 
C-|—|- application and can be found in [4]. The program is able to find all invariants of the 
MSSM (derived in section 3) within seconds. Additionally, two more invariants have been 

®Such a restriction could also be implemented by for example including negative powers. 
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found for the dMSSM: 


11 7 = 5 '27501^-31965^25920 (y„^y„^y„3)-3859 {y ^^51 (4 

1 18 = 309Mi + 4059M2 - 6336M3 - 693 Tr - 495 Tr {adVcT^) 

- 242 Tr (tteye"^) + (4.18b) 

For more details on the implementation of the above method, as well as a further extension, 
see [5]. 


4.3 Higher loop orders 

The computer algebraic method can easily be extended to higher loop orders. To this 
end, let us consider the general form of a /3-function for Xi in terms of its different loop 
contributions: 

/3(x,) = (x.) + ^/3(2)(x,)+ ... (4.19) 

Two-loop invariants can therefore be found by considering candidate invariants of the form 

/ = + (4.20) 

The derivative of I with respect to t reads 


— = + — 
dt ^ IGvr^ 






(IGvr 


1_ r(2) 

2 i 2 -'2 > 


(4.21) 


where I) 2 is the contribution to the derivative that involves the j-th loop order /3-function 
Moreover, the terms have been grouped by equal powers of the factor l/(167r^). 
The requirement for RG invariance now reads: 

/W = 0, if) + if) = 0, (4.22) 


which can be reduced to a linear system of equations, equivalent to the method for one-loop 
invariants. 

The computer algebraic method has been applied to the one and two-loop /3-functions 
of the MSSM, the dMSSM, and the pMSSM. All one-loop invariants found are consistent 
with the ones constructed in section 3. While the one-loop /3-functions are quite well-known, 
the two-loop ones are not. They were taken from [6] and thoroughly checked against the 
results of [7-9]. For the MSSM, we have found the following two-loop invariant: 

Ji = ^ (^Mi + 209M2 - 88M3 + . (4.23) 
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Note that, apparently, only the one-loop invariant M 2 /g^ has a two-loop continuation in 
the MSSM. For the dMSSM, we have found two additional two-loop invariants: 


363Mi 1 

•^2 = -P2 ^ 




g'^ 1671-2 

- 1243Tr (adi/d“^) 
IIM 3 1 


894Mi 6732M2 - I 6 IO 4 M 3 - fill Tr 

590761 




gj 167r2 


336 


66 M 3 -h llTr -h llTr (adi/d - — 


and for the pMSSM we have found 


(4.24a) 

(4.24b) 


J 4 = -( 2869Mi 1485M2 - 13640M3 3762^7 -h 3498^^ 

g'^ lOTT^ \ 

_ 15186\ 

^ ) ’ 

J5 = -( 227 Ml + 386IM2 - 3586M3 - 198^7 - 330^^ ) . 

5s 167r2 V /i y 


(4.25a) 

(4.25b) 


Due to simplifications in the trilinear sector with respect to the MSSM, the dMSSM and 
pMSSM have two-loop continuations of the one-loop invariants Mi/g''^ and Mj,lg^ as well. 


5 Conclusion and outlook 

We have developed two novel, efficient methods for finding RG invariants. The more 
theoretically inclined approach links the existence of one-loop invariants to symmetries of 
the underlying theory. For any given supersymmetric theory that has the same structure 
as the MSSM, the number of RG invariants that involve scalar masses is equal to the 
number of U(l) symmetries of its family sector. The computer algebraic method is able to 
find invariants at higher loop orders and is applicable to any set of RG equations. Both 
methods have been applied to the /3-functions of for example the unconstrained MSSM and 
the pMSSM. For the MSSM, three new invariants at one-loop order, and one new invariant 
at two-loop order have been found. For the pMSSM we have found three new invariants 
at two-loop order. 

A next step in the development of using RG invariants as probes of high-scale physics, 
could be the construction of new sum rules for various supersymmetry breaking models. 
This has been done in [3] for the previously known invariants in the pMSSM, and could 
now be extended to the (d)MSSM. In case supersymmetry is found, it is unlikely that the 
entire spectrum will quickly be measured. Even our present knowledge of the Higgs mass 
and other supersymmetry-sensitive data only serves to constrain a number of parameters, 
depending on the supersymmetry breaking scenario. However, one of the advantages of 
using RG invariants to probe high-scale physics is the fact that not all parameters of the 
theory need to be included, simply because the invariants typically contain only a subset 
of all the parameters, and some parameters do not enter at all. 

Furthermore, since most invariants are directly linked to symmetries of the underlying 
theory, there is a certain amount of freedom to choose what parameters to include in the 
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invariants. Thus, one might be able to exclude certain parameters from part of the analysis 
(i.e. particular sum rules of invariants) by cleverly picking linear combinations of quantum 
numbers such that those parameters do not appear in the relevant invariants. Of course, it 
could turn out that physics beyond the Standard Model matches a different (perhaps non- 
supersymmetric) effective field theory rather than one that is discussed in this paper. Since 
the RG invariants method is completely general, it could still be used to probe high-scale 
physics models once the /^-functions of the appropriate effective field theory are known. 


A The MSSM 


After we consider a general supersymmetric Lagrangian, we will provide the held content 
and Lagrangian for the MSSM, its U(l) symmetries, and some one-loop /3-functions. 


A.l A general supersymmetric theory 

In supersymmetric theories the helds are grouped into supermultiplets. A chiral supermul- 
tiplet (labeled by i, j) consists of a complex scalar held (j) and a left-handed Weyl spinor ijj. 
Each gauge group (labeled by v) with corresponding inhnitesimal generators r“, structure 
constant and gauge coupling gives rise to a gauge supermultiplet. A gauge super- 
multiplet consists of real gauge boson helds and Weyl spinors A“. The gauge invariant 
superpotential is denoted by W, and we dehne 



tRb' = 


6‘^W 

64>i6(j)j' 


(A.l) 


Using this notation, a general supersymmetric Lagrangian then has the following form [10]: 
^SUSY = - 

- v /2 • K + h.c.) - |iuf - ^ + h.c.) 

(A.2) 

where = (l, —cr^, —cr^, —cr^) involves the unit matrix I and the Pauli matrices (t“ . The 
indices of Weyl spinors are raised and lowered by the antisymmetric symbol with non¬ 
zero components = —ei 2 = 621 = 1. By convention, spinor indices are always 

contracted diagonally downwards for left-handed Weyl spinors. The dot in eq. (A.2) is 
used to denote the contraction of two Weyl spinors that yields a spin-0 singlet.The held 
strength tensors are dehned by 


= f) 8 4°- — n AC 

VflU — yvJv 


(A.3) 


^To avoid clutter, though, this dot is implicit in the remainder of this appendix. 
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and the gauge covariant derivatives read 


{D^X^r = AS, 

(A.4a) 

D^cpi = {df, -b ig^Al^T^) cpi, 

(A.4b) 

D^ipi = {df, -b igvA^^rpp) 'ipi. 

(A.4c) 


In the remainder of this appendix, the gauge couplings belonging to the gauge groups 
U(l), SU(2), and SU(3) are denoted by g', g, and gs respectively. The completely anti¬ 
symmetric structure constants are respectively given by 0, , and (the primes are 

used to distinguish SU(2) from SU(3) indices). The generators of SU(2) and SU(3) are 
proportional to the Pauli matrices (t“ and the Gell-Mann matrices A“ respectively. 

A.2 The MSSM Lagrangian 

The field content of the unbroken MSSM, including the corresponding gauge group rep¬ 
resentations, is given by tables 5 and 6. To distinguish the notation for Standard Model 
fields from their supersymmetric partners, the latter receive a tilde. In these tables the 
right-handed Weyl spinors have been conjugated to bring them in a left-handed form: for 
a Weyl spinor ijj we define 


Name 

Spin 0 

Spin 1/2 

SU(3)c X SU(2 )l X U(1)y 

sleptons, leptons 

Ll = {^L Gl) 

Ll = {i^L ei) 

(1,2,-i) 

(1,1,1) 


Ql = {UL di) 

Ql = {uL di) 

(3,2,1) 

squarks, quarks 

~5k 

Ur 

Ur 

(3,1,-i) 


^*R 

^R 

(3,l,i) 

Higgs, Higgsinos 

= (H* H») 
Hi = ( Hi HR 

= (tr+ #2) 
Hi = (H'i HR 

(l,2,i) 

(1,2,-1) 


Table 5. Chiral supermultiplets of the unbroken MSSM and the corresponding gauge group rep¬ 
resentations. 


Name 

Spin 1/2 

Spin 1 

SU(3)c X SU(2)l X U(1)y 

gluinos, gluons 

r 


(8,1,0) 

winos, W bosons 

IT“' 


(1,3,0) 

bino, B boson 

B 

B 

(1,1,0) 


Table 6. Gauge supermultiplets of the unbroken MSSM and the corresponding gauge group rep¬ 
resentations. 
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Including the phenomenologically motivated requirement of i?-parity conservation, the 
MSSM superpotential is given by [10]: 

H^MSSM = - e^jiye{LL)°‘{Hd)a + “ dj^VcliQ lT {Hd)a 

+ fi{Hur{Hd)a. (A.5) 

Note that in this expression all color and family indices have been suppressed. For conve¬ 
nience, we do show explicitly the weak isospin doublet indices that are raised and lowered 
by €^13 (like the Weyl spinor indices). All parameters of the MSSM are defined in table 7 
at the end of this subsection. 

The full Lagrangian of the MSSM (including soft supersymmetry breaking terms), split 
up in parts, is given by [10]:^^ 


• Kinetic terms for the gauge supermultiplets and gauge interactions: 


MSSM ^ 


- (A.6) 

where the gauge covariant derivatives for the gauginos read 

(D^gT = - gsr’^^G’lW, (A.7a) 

{D^wf = (A.7b) 

The field strength tensors are given by 

G% = d^Gl - d,Gl - g^r^^GlGl, (A.8a) 

- d,W-' - ge-'^'^'W^JW^\ (A.8b) 

Bf,u = dfj,By - duB^. (A.8c) 




• Kinetic terms for the chiral supermultiplets and gauge interactions: 




MSSM 


{D^^Ll)^D^Ll + iL{a^^D^LL 
+ {D>^rn)^D^ek + iep^D^e% 

+ {Df^QL)^D^QL + iQ{a^^D^QL 
+ D -k D 

+ {D^d},)^D^d*j, + idp^D,d% 

+ + iHla>^D^Hu 

+ {D^Hd)^D^Hd + iHyD^Hd, 


(A.9) 


11 


As gauge fixing terms are irrelevant for our analyses, we simply omit those here. 
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where the gauge covariant derivatives for the leptons, quarks and Higgs doublets are 
given by 


D^Ll = ( 5 ^ + \igW^^a^' - \igB^ U, (A.lOa) 

= (^5^ + ig'B^ e% (A.lOb) 

D^Ql = ( 5 ^ + \ig.Gl\^ + \igW-^<y^' + Ql, (A.lOc) 

- lig'B^^ u% (A.lOd) 

- \wsGl\-* + \i9'B^ d% (A.lOe) 

D^Hu = ( 5 ^ + \igW^^a<^' + \igB^ (A.lOf) 

+ ligW^'a^^' - ligB^^ H^. (A.lOg) 


As superpartners have equal quantum numbers, the covariant derivatives for the 
sleptons, squarks, and Higgsinos are exactly the same. 

• Chiral supermultiplets coupled to gauginos: 

-^/2 [g.4T:^. ■ A“ + h.c.)]^^^^ = - ^/2 ^^gL[a-'Lr.W^' - \g'L{LLB 

+ g%e<RB + ]^g,QlX^QLT 

+ \9QWQlW‘^' + \gQ{QLB 
~ — -g'u^u^B 

~ 29sdRy*<fRg°' + -g'd-Rd^B 
+ ^gHla^'HuW^' + ^g'HlHuB 
+ ^gHy'HdW^^' - \g'HyB + h.c.) . 

(A.ll) 


- 22 ~ 



• Scalar interactions coming from the “F-fields”: 

- I^TJmssm ^ “ L\yleRe\yJ.L - e^RyeyleRH^Hd 

+ {H*df{V^)l3ylye{LLT{Hd)a “ Q^vIuR^^V^Ql 

- u^^yuyiuRHiH^ + {Kf{Qi)pyiyu{QLr{Hu)o. 

- QWd^RdRVdQL - dlydy^dRE^Hd 

+ {H*df&L)pydydi.QLT{Hd)a - (Q{y\dRe^jiyJ^L 
RRyuV(I^ rH d^u d' ^Ryedd^L/R /i uj^yuHdQL 

- y*dlydHiQL + h.c.) - \yi\‘^HlHu - (A.12) 


Yukawa couplings: 


1 


-- [W^^i • + h.c.) 


J MSSM 


e^Rye{LLr{Hd)a + e^Rye{LLT{Hd)c. 

+ e^^ye{LLT{Hd)a - U^Ryu{QLT{Hu)a 
- U^Ryu{QLT{Hu)a - u'Ryu{QLT{Hu)a 
+ d}Ryd{QLT{Hd)c + d^Ryd{QLT{Hd)a 
+ dlyd{QLT{Hd)a - yi{HuT{Hd)c. + h.c. (A.13) 


Four-scalar interactions coming from the “D-fields”: 


1 


V \r\ddv(t>i 


J MSSM 


(qIa“Ql - u^r\^ur - 4A “ dH )' 

^5^ {lia^'lL + Q4“'Ql + + Hy Hd] 


- y ( - e\eR - y\QL + yRUR 


- \dldR - IhIHu + \HlHd 


(A.14) 


Soft supersymmetry breaking terms: 
1 


^MSSM = - § (+ M 2 lY“'' 1F“' + MiB^B + h.c. 




a'T^a' 

■t ^ (J^\0‘(TT\ _ ~t „ (n^\ot(TT\ iJt, 


+ 


e\ae{LLT{Hd)a - u\au{QLT{Hu)a + d]^ad{QLT{Hd)a + h.c. 

- L\rri~LL - 

- miyHu - ml^Hd - [b{H^r{Hd)a + h.c.] . (A.15) 

All parameters of the MSSM are listed in table 7. Many degrees of freedom are 
unphysical though, as they can be absorbed by clever field redefinitions.^^ 


^^The total number of independent parameters in the MSSM is 123 (not including the strong CP violating 
angle). 


- 23 - 











Name 

Physics description 

Math description 

# 

gs,g,g' 

gauge couplings 

real numbers 

3 

Ve-! Vu"} Vd 

Yukawa coupling matrices 

complex 3x3 matrices 

54 

Mi,M2,M3 

gaugino masses 

complex numbers 

6 


trilinear coupling matrices 

complex 3x3 matrices 

54 

9 9 9 9 9 

mi, m|, mi, mi, mi 

sfermion mass matrices 

Hermitian 3x3 matrices 

45 


Higgs masses 

real numbers 

2 


Higgs mixing parameters 

complex numbers 

4 

Total: 

168 


Table 7. All parameters of the MSSM including soft supersymmetry breaking terms. The param¬ 
eters that are listed above the horizontal dashed line also occur in the Standard Model. 


A.3 U(l) symmetries 

The MSSM Lagrangian including soft supersymmetry breaking terms has three U(l) sym¬ 
metries: weak hypercharge (Y), baryon number (B), and lepton number (L). The last two 
symmetries are present as a consequence of imposing i?-parity conservation. If we ignore 
the soft supersymmetry breaking terms and set the supersymmetry preserving Higgs mix¬ 
ing parameter to zero, then, as it turns out, there are two more U(l) symmetries in the 
MSSM. The corresponding groups, that we will call U(l)yi and U(1 )b, are equivalent to 
the Peccei-Quinn (P) and R symmetries (see [11, 12]). Let (p denote any field contained in 
the MSSM, then the transformations under U(l)y 4 and U(1 )b are respectively defined as 

P —^ (p —^ (A. 16) 

The charges Qa, Qb depend on the fields and uja,ojb are free parameters. Tables 8 and 9 
list the values of Qa, Qb for the fields that make up the chiral and gauge supermultiplets 
respectively. 

Now let us consider the MSSM in full (including the soft supersymmetry breaking 
terms), and let us apply the U(l)y 4 and U(1 )b transformations to all fields. If we would 
require invariance of the MSSM Lagrangian under these transformations, then we would 
need to simultaneously redefine the Higgs mixing parameters as follows:^^ 

^ ^ ^ (A.17) 

and the gaugino masses and trilinear coupling matrices as 

Mfc ^ (A.18) 

the transformations under U(1 )a and U(1 )b are parametrized by two independent parameters 
ujAyijJB, two objects in eqs. (A.17) and (A.18) can each get one of their phases removed by fixing uja,^b 
in a clever way. Conventionally, the soft breaking parameters b and Ms are made real in this way. 


- 24 - 











Spin 0 

Qa 

Qb 

Spin 1/2 

Qa 

Qb 

Ll 

0 

1 

Ll 

1 

0 


0 

1 

^R 

1 

0 

Ql 

0 

1 

Ql 

1 

0 

^*R 

0 

1 


1 

0 

d* 

0 

1 

^R 

1 

0 

Hu 

-2 

0 

Hu 

-1 

-1 

Hd 

-2 

0 

Hd 

-1 

-1 


Table 8. The charges Qa, Qb for the fields that make up the chiral supermultiplets. 


Spin 1/2 

Qa 

Qb 

Spin 1 

Qa 

Qb 

r 

-1 

1 


0 

0 


-1 

1 


0 

0 

B 

-1 

1 

B 

0 

0 


Table 9. The charges Qa, Qb for the fields that make up the gauge supermultiplets. 


where k = 1,2,3 and '0 = e,u,d. From eq. (A.18) we infer that the family sector of 
the MSSM (which does not include /i and b), is not invariant under U(l)yi and U(1 )b 
separately, but only under the combination U(1 )a+s-^‘^ This means that on top of Y, B, L, 
there is a fourth quantum number X = Qa + Qb that corresponds to a U(l) symmetry of 
the MSSM family sector. 

A.4 One-loop /3-functions in the MSSM 

The one-loop /3-functions for the gauge couplings are given by 

f3{g') = ll5''^ I3ig) = g^, /3(5s) = -^gt (A.19) 

and those for the squares of the gaugino masses read [7] 

(3 (|Mi|2) = 445'2|Mi|2, /3 {\M2f) = (3 {mf) = -l2g^,\Msf. (A.20) 

^^One may argue that p is contained in the family sector of the theory; this is a matter of taste. However, 
since /r is a supersymmetry preserving parameter, it turns out to be completely irrelevant for our analyses 
in section 3. Hence, we can safely ignore this parameter and simply define the MSSM family sector to not 
include p. 


- 25 - 





















The one-loop /3-functions for the soft scalar masses are given by 

= - 8y|^5'2|Mip - - y + 2m‘jj^ylyu + 2m\^ylyd 

+ 'I^^ylvu + m'^Qylyd + ylz/um| + ylyd.m'^Q + 2ylmly^ 

+ 2y\m‘^~yd + 2a]^a„ -b 2a\ad + 2YQ^g'‘^S, (A.21a) 

/3(m|) = - 8Y~^g'^\Mi\‘^ - dg^\M2\‘^ + 2m|^^r/|ye + rri~ylye + ylyem^ 

+ 2ylm~ye + 2a^ae -b 2Yj^J'^S, (A.21b) 

32 

/3("ii) = - + Am^j^yuyi + ‘^^ly^yl + ^y^^yl^l 

+ Ayum^yl + 4a„a]) -b 2Yu*J‘^S, (A.21c) 

= - 85^1 - ‘^gl\M^\^ + Amjj^ydyl + + 2ydl/d«^| 

R O 

-b Ayd-m^yl -b AadO^ -b 2Y^,j)'‘^S, (A.21d) 

/3(m|) = - 8Y~,^g‘^\Mi\^ + Am]j^yeyl + ‘^^^y^yl + 2yeylml + Ayem^yl 

+ Aaf,a^ + 2Ygt^g'‘^S, (A.21e) 

+ QYT[rn]j^ylyu + m^^ylyu + ylrnly^ 

+ 4 a„) + 2YhJ^S, (A.21f) 


- 65^1^2^ -b 2Tr(^3m|^^4l/d + 3m|4yd + mjj^ylye 

+ m^ylye + 2>ylrin?~yd + ylmlye + -b a^ae) + 2YH^g'‘^S. (A.21g) 

The quantity S that appears in the /3-functions above arises from tadpole diagrams 
and is defined as 

j) = ^ ("^1 “ "^1 “ + m| + m|) + . (A.22) 

^ 4> ^ 

The /3-function for S follows directly from the /3-functions for the soft scalar masses and is 
given by 

/3(5) = 2 ^ = 2245 . (A.23) 

0 

The /3-functions for the soft scalar masses that are given above are in correspondence 
with the ones in [7]. The advantage of our result though, is that these /3-functions are 
expressed in terms of the weak hypercharges of the fields, which is useful for our symmetry 
analyses. 

B The dMSSM 

The dMSSM is the “flavor-diagonal” version of the MSSM. The simplifications with respect 
to the MSSM are as follows: 
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The Hermitian sfermion mass matrices m ~ with f = L,e, Q, u, d are taken diagonal, 


i.e. 




/ m% 0 0 \ 

h 

0 m~ 0 
/2 

0 0 m~ 

73 / 


V 


(B.l) 


The Yukawa and trilinear coupling matrices and with ^ = e,u,d are taken 
real and diagonal, i.e. 


(B.2) 




0 

»1 


0 

0 ^ 

Vip — 

0 

y'i/>2 

0 1 ? — 

0 

0.02 

0 


1 0 

0 

y-ips / 

1 0 

0 

003 J 


All parameters of the dMSSM are listed in table 10. 


Name 

Physics description 

Math description 

# 

9s,g,g' 

gauge couplings 

real numbers 

3 

Vei Vu: yd 

Yukawa coupling matrices 

real, diagonal 3x3 matrices 

9 

Ml, M 2 , M 3 

gaugino masses 

complex numbers 

6 


trilinear coupling matrices 

real, diagonal 3x3 matrices 

9 


sfermion mass matrices 

real, diagonal 3x3 matrices 

15 


Higgs masses 

real numbers 

2 

y,b 

Higgs mixing parameters 

complex numbers 

4 

Total: 

48 


Table 10. All parameters of the dMSSM including soft supersymmetry breaking terms. The 
parameters that are listed above the horizontal dashed line also occur in the Standard Model. 


C The pMSSM 


The pMSSM is a heavily simplified version of the MSSM. Some parameters in the MSSM 
give rise to processes that seem improbable from a phenomenological point of view, such 
as flavor-changing neutral currents and CP violation beyond experimental bounds. To 
suppress these possibilities, one usually imposes several constraints. For the pMSSM they 
are as follows: 


The Hermitian sfermion mass matrices m ~ with f = L,e, Q, u, d are taken diagonal, 
and the first and second generation masses are assumed to be degenerate, i.e. 



/ m% 0 0 \ 

ji 

0 m~ 0 


(C.l) 
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• The first and second generation Yukawa couplings are neglected and the third com¬ 
ponents are taken real: 



/o 0 0 \ 


/o 0 0 \ 


/o 0 0 \ 

Ve = 

0 0 0 

) Vu = 

0 0 0 

, Vd = 

0 0 0 


0 Vr J 


\0 0yt j 


Q ybj 


• The trilinear coupling matrices are taken proportional to the corresponding Yukawa 
coupling matrices, which implies 

CLt — -^Tyri CLt — CLb ~ -^bUb- (C'^) 

• The gaugino masses and Higgs mixing parameters are taken real. 

All parameters of the pMSSM are listed in table 11. 


Name 

Physics description 

Math description 

# 

9s,g,9' 

gauge couplings 

real numbers 

3 

yT,yt,yb 

Yukawa couplings 

real numbers 

3 

Ml, M 2 , M 3 

gaugino masses 

real numbers 

3 


trilinear couplings 

real numbers 

3 

9 9 9 9 9 

sfermion mass matrices 

real, diagonal 3x3 matrices 

10 


Higgs masses 

real numbers 

2 

y,b 

Higgs mixing parameters 

real numbers 

2 

Total: 

26 


Table 11. All parameters of the pMSSM including soft supersymmetry breaking terms. The 
parameters that are listed above the horizontal dashed line also occur in the Standard Model. 


D Summary of results 

This appendix lists all RG invariants that have been found for the MSSM, the dMSSM, 
and the pMSSM. The one-loop invariants have been found by both using the approach 
involving symmetries and by applying the computer algebraic techniques. The two-loop 
invariants, however, have only been found by the latter. 

D.l The MSSM 

Tables 12 and 13 contain all one and two-loop RG invariants respectively that have been 
found for the (unconstrained) MSSM. 
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# 

Q 

RG invariant 

1,2,3 


^1 ^2 ^3 

~9F 

4,5 


j_ 11 jy 1 n 

gl2 g2 ) gl2 “T ^2 

6 

Y 

1 

IP 

Tr (m ~ — m~ — 2m~ -|- m ~ -|- m~ ) -|- rtiir — 111 % 

\ Q L ^ d ^dj 

7 

- L)-8Y 

Tr ( 14m ~ — 14m~ -|- 5m~ — 19m ~-|- 3m?^ — 8m%, -|- 8m^ 

\ Q L ^ d Hu ' Hi 

8 

3B + L 

Tr (^6m~ -|- 2m|. — 3m~ — 3m~— m|^ — Mip -|- 36 M 2 P 

9 

X 

Tr 1 

6mT -|- 2m~ -|- 3m~ -|- 3m~ -|- m~ | — 4m^ — 4m^ 

, Q L u ' d ej Hu Hi 

+ ^|Mi|2 + 24|M2|2-32|M3|2 


Table 12. One-loop RG invariants in the MSSM. The second column, if applicable, lists the 
quantum number Q that corresponds to the invariant. 


# 

RG invariant 

1 

11 M 2 1 

167r2 

'Ml + 2 mM2 - 88M3 -b 



Table 13. Two-loop RG invariant in the MSSM. 


D.2 The dMSSM 

Tables 14 and 15 contain all one and two-loop RG invariants respectively that have been 
found for the dMSSM. 


D.3 The pMSSM 

Tables 16 and 17 contain all RG invariants that have been found for the pMSSM. These 
invariants are also listed in [1-3], but often as different linear combinations. 


- 29 - 




















Table 14. One-loop RG invariants in the dMSSM. The second column, if applicable, lists the 
quantum number Q that corresponds to the invariant. 


I 









Table 16. One-loop RG invariants in the pMSSM. The second column, if applicable, lists the 
quantum number Q that corresponds to the invariant. 




# 

RG invariant 

1 

^ (^2869Mi -h 1485M2 13640M3 -h 3762A -b 3498A^ 

2 

^ (Mi + 209M2 - 88 M 3 + 

3 

16 ^ 2 ( 227 M 1 -b 386 IM 2 3586 M 3 198A 330A^ -b 


Table 17. Two-loop RG invariants in the pMSSM. 
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